Self- Affinity in the Gradient Percolation Problem 
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We study the scaling properties of the solid-on-solid front of the infinite cluster in two-dimensional 
gradient percolation. We show that such an object is self afBne with a Hurst exponent equal to 2/3 
up to a cutoff-length ~ 5 where g is the gradient. Beyond this length scale, the front position 
has the character of uncorrelated noise. Importantly, the self-afflne behavior is robust even after 
removing local jumps of the front. The previously observed multi afflnity, is due to the dominance 
of overhangs at small distances in the structure function. This is a crossover effect. 
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Rough surfaces showing non-trivial scaling properties 
have been extensively studied theoretically, numerically 
and experimentally over the last couple of decades. Ex- 
amples of such surfaces are those appearing during brittle 
fracture which were first characterized as being frac- 
tal 01 but it was then realized that the concept of self 
affinity was more appropriate The question of self 
affinity versus fractality has also been the focus of in- 
tense research on invasion fronts in porous media and on 
the dynamics of magnetic domain walls 0|. It was re- 
cently reported that the displacement fronts in self-affine 
fractures are self affine 0. More recently, a possible ex- 
planation for the observed self affinity of fracture surfaces 
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FIG. 1: Top-side and bottom-side SOS fronts based on the 
perimeter of the cluster connected to the p = 1 edge (shown 
in the insert as dots). We also show the filtered jo{i) front 
(see Eq. ®). 
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has been proposed and hinges on a clear understandingof 
the distinction between fractality and self affinity 0, □] ■ 
It has also been suggested that brittle fracture surfaces 
are multi affine rather than simply self affine ■ Whether 
this is so remains an open question j^. 

It is the aim of this Letter to study the question of 
fractality, self affinity and multi affinity of a front in a 
system which is simple enough to be tractable, namely 
that of the gradient percolation There are already 
in the literature studies of this system in the present con- 
text. Furuberg et al. m study the jumps in the position 
of the solid-on-solid (SOS) front of the infinite cluster, 
whereas Asikainen et al. |0| conclude that this front 
is multi affine. We will in this Letter show that up to a 
given scale, the SOS front is self affine with a well-defined 
Hurst exponent, whereas on larger scales its position be- 
comes uncorrelated. The self affinity is not caused by the 
jumps in the position of the front due to overhangs, but 
related to its fractal structure. The multi affinity seen by 
Asikainen et al. has its origin in the overhangs resulting 
from the definition of the SOS fronts and shows up in the 
structure function on small scales. 

In gradient percolation, a spatial gradient in the oc- 
cupation probability p is introduced. A cartesian coor- 
dinate system is oriented with respect to the finite 
lattice of size Li x Lj (assuming for the rest of this pa- 
per that the lattice is two dimensional), so that the i 
axis runs perpendicular to the gradient (i.e. along the 
lower edge) and the j axis along the gradient (i.e. the 
left edge). The gradient is introduced in the j direction 
so that p{j) = gj, where the gradient g = How- 
ever, the cluster connected to the lower edge will reach 
some average value, j ~ jg, with an associated occupa- 
tion probability pg — gjg . The region around jg is critical 
and has a width ^, spanning between j± = jg±^/2, where 
^ is the correlation length associated with the critical re- 
gion in the direction of the gradient. Defining p± = gj± 
and setting £, = \p± - Pgl^" = \g{j± ~ jg)]'" wlwre v is 
the correlation length exponent, Sapoval et al. [l^ found 
that C ~ = <7-4/7. 
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FIG. 2: Data collapse of the averaged wavelet coefficients for 



the bottom side front based on lattice size L, 



64 to 8192, 



while Li = 2048. We have that g = 1/Lj. 
has a slope of + 1/2, see Eq. 



The straight line 



The infinite cluster has a fractal structure "with an up- 
per cutoff in length scale set by the -width of the critical 
region, ^. We now focus on the front of this infinite clus- 
ter and define precisely "what we mean by this front 
in the gradient percolation problem. Our starting point 
is the perimeter of the cluster of occupied sites that is 
attached to the p = 1 edge of the lattice. Since this 
perimeter contains overhangs and therefore is multival- 
ued "when interpreted as a function "we use the SOS 
method to extract a single-valued function for its posi- 
tion, see Fig. For each i, we use either the j value 
that is closest to the p = edge (top side) or the j value 
"which is closest to the p = 1 side (bottom side) or the 
average over all the j values attached to a given i value 
(average front). 

A trace is statistically self affine if the probability 
density, 7r(z,j), for it to have a value j at i, given that 
j = at I = 0, has the invariance 



A'^7r(Ai, A'^j) =7r(i,j) , 



(1) 



"where C is the Hurst exponent. This invariance must be 
caused by spatial correlations in j along the i-axis. We 
note that a Levy flight, which is an uncorrelated random 
walk whose step size h is drawn from a power law distri- 
bution N{h) ~ h-^-'^, will satisfy Eq. ^ with an appar- 
ent Hurst exponent C = However, in this example, 
satisfying Eq. is due to the step size distribution and 
not to spatial correlations [l5| . 

We have used the Average Wavelet Coefficient (AWC) 
method fvA IT^ to analyse the structure of the SOS 
fronts. The AWC method consists of wavelet transform- 
ing and averaging the wavelet coefficients w{b, a) at 
each length scale a over position 6, W{a) = {w{b,a))b. 
If j{i) is self affine, the averaged wavelet coefficients will 
scale as 



W{a) 
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C+l/2 
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FIG. 3: Data collapse of the averaged wavelet coefficients for 
the smoothed io(«) based on the bottom side front. The lattice 
sizes and gradients are as in Fig. |21 The long-dashed line has 
a slope of 7/6, see Eq. whereas the dotted line has a slope 
of 1=1/2+1/2, consistent with uncorrelated random walks. 



We show in Fig.|21 the averaged wavelet coefficients based 
on the Daubechies-4 wavelets for the bottom-side fronts. 
The plots for the top and average fronts are comparable. 
The data are based on averages over 2201 samples for Lj 
in the range 64 to 2048 and 200 samples for Lj = 4096 
and 8192. L^ was set to 2048 for ah the different Lj. The 
gradient g was set to l/ij. There is a clear crossover 
between two regimes in these plots. At smaller length 
scales, one does indeed find the behavior of Eq. ^ in- 
dicating self affinity. On larger scales, the slope of the 
log-log plots are zero indicating ^ ~ —1/2, which corre- 
sponds to uncorrelated or white noise |19| . Furthermore, 
we observe excellent data collapse when W is scaled by 
g^^ and the length scale, a, is scaled by We will 

show below that 



and 
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where D^. = 4/3 is the fractal dimension of the external 
perimeter of the front plj . 

The main goal of this Letter is to derive Eq. Q and 
thus demonstrate that ^ is a proper Hurst exponent and 
i{i) a sclf-affine function. To this end, we need to demon- 
strate two things: First, j{i) satisfies the scaling relation 
and, secondjthat this is not due to a power law tail in 
the step size distribution. We note that since the aver- 
age wavelet coefficients obey Eq. 0, j{i) automatically 
satisfies Eq. (Q. Therefore, we now need only to identify 
the mechanism behind this scaling. 
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In order to derive Eq. Q, we start by noting that the 
distribution of distances m between crossing points be- 
tween a planar fractal curve with dimension D^. — e.g., 
the percolation perimeter — and a straight line follows 
the power law 7r(m) m~^= |20l |. Introducing a gradi- 
ent in the j direction and placing the straight line in the 
critical region interval, [j-, j+J, and parallel to the i axis, 
the distribution of crossing point distances m remains the 
same. A self-afhne curve characterized by a Hurst expo- 
nent C, leads to a distribution of crossing point distances 
given by 7r(m) ~ m~(^~^) pol |. By comparing this ex- 
pression to 7r(m) ~ m"^", Eq. ||2Jl immediately follows. 
However, we still need to show that is indeed self 
affine, in other words the scaling relation Eq. is not 
caused by jumps. 

First, we turn to deriving Eqs. and jSJ. The cor- 
relation length in the direction of the gradient, the j di- 
rection, is ^ ~ g-'^/(i-+<^) _ Since the perimeter is locally 
isotropic, this is also the correlation length in the i direc- 
tion. The crossover length scale from self affinity to un- 
correlated noise is the correlation length ^. Hence, rescal- 
ing a — >• a/^ ~ a/ g~'^^^^~^'^^ gives data collapse along this 
axis which demonstrates Eq. Q). Likewise, the crossover 
length scale in the j direction is ^. This implies that the 
normalized wavelet coefficient at this scale, W{S,)/^^^^ is 
equal to Hence, W{^) ~ ^3/2 _ ^-(3/2)0 _ ^-/3^ 

P = (3/2)a = 6/7, as stated in Eq. ®. 

In order to show that is a self-affine function, we 
need to demonstrate that C is not caused by the step size 
distribution. To this end, wc define the following trans- 
formation of the function jk{i) where we factorize 
the function in such away that we can distinguish the 
respective roles of persistency and step sizes, 

i 

jk{i) = ^ sign[j(m+l)-j(m)] \j{m+l)-j{m)\'' , (6) 

where \j(m+l)—j(rn)\ = h(rn) is the step size at position 
m. We have in particular that = It shown 
in [T]| , that h is distributed according to 

N{h,g)^h-'''-'f{hg^) , (7) 

where = 4/3 and /(z) approaches a constant as z <C 1 
and falls off faster than any power law as z — > cx). The 
step size distribution comes from the appearance of over- 
hangs in the perimeter. An overhang is defined as the 
jump made by the front from one position along the i axis 
to the next due to a backwards turn [TlL[T^[T4| . In order 
to confirm that the overhangs do not generate the Hurst 
exponent ( = 2/3, wc analyse the filtered front jo{i), de- 
fined in Eq. 1^. With fc = 0, we eliminate the overhangs 
all together |22|. Fig.|31shows the data collapse based on 
jo{i) corresponding to the bottom side shown in Fig. 
|5] The scaling along the i axis is unchanged as no change 
in the system has been made in that direction. However, 
since all step sizes have been reset to unity in the trans- 
formation — > jo(*)j the rescaling in the j direction is 



no longer controlled by jc- In order to regain data col- 
lapse for different g = l/Lj, we need to rescale the lattice 
units in this direction by the Hurst exponent, C = 2/3. 
The straight line matching the small-a region of the fig- 
ure has a slope 2/3 + 1/2, while the straight line matching 
the large-a portion has a slope of 1/2 + 1/2 correspond- 
ing to an uncorrclatcd random walk. This shows that, 
for small scales, the C = 2/3 is indeed a Hurst exponent. 
On the other hand, for longer length scales, we expect 
random walk behavior since white noise gives precisely 
the exponent 1/2 in the transformation —^ joii)- 

In order to analyse the multi affinity that has been 
reported in this problem |l2l |. we construct the struc- 
ture function Ck{n,g) = {\j{m + n) ~ j(m)|'''). Multi 
affinity occurs when Ck{n, g)^^^ does not scale with a 
single fc-indcpendent exponent with respect to n. Using 
the overhang distribution (U)), we find Cfe(l,g) ~ g'^''^\ 
where s(fc) = ijmi[Q , a{D ^ - k)] = min[0, (16/21 -4A;/7)]. 
The self-affine character of cannot be visible in the 
structure function for n = 1 but will appear only grad- 
ually as n is increased. We may therefore analyse the 
structure function based solely on the Levy character in- 
duced by the overhangs in the small-n limit. We will 
call this the Levy regime, whereas for larger n where the 
self affinity dominates, we will refer to as the self-affine 
regime. The scaling with respect to g for Cfe(l, g) persists 
for n > 1 in the Levy regime since j{i-\-n) follows a 

Levy distribution whose power law tail does not change 
with increasing n. Hence, we expect Ck{n,g) ^ in 
this regime. In order to derive its dependence on n in 
the Levy regime, we note that the distribution of dis- 
tances I between overhangs follows the same power law 
as the overhangs themselves. This can be seen as follows. 
When there is a gradient present in the j direction, the 
length of the perimeter scales as L^", when the gradi- 
ent is kept fixed. Making a cut through the perimeter 
with a straight line parallel to the i axis, the crossing 
points of the perimeter with the line form a fractal set 
with dimension — 1. Hence, there are, in a given in- 
terval I, Ni ^ l^""^ overhangs |23|. These overhangs 
give rise to an effective Hurst exponent l/D^ = 3/4 
on the fractal set, seen e.g. in the width of the trace, 
Aj ~ Ni^^' ~ i{Dc-i)/Dc ^ Since the overhangs form a 
fractal set, we will need Ni, ~ ;-(^<!-i) boxes of size I to 
cover it. Due to the averaging over position i, there will 
be yet another factor see [23 • We may now assemble 
these pieces to form the scaling of the structure function 
in the Levy regime, Ck{l,g) Nj^Nbl ~ I'^'f^ where 

.r \ 1 1 2-i:>e 1 2 

Therefore, in the Levy regime, i.e. for small n, there is 
multi affinity. A similar analysis in the self-affine regime, 
i.e. at larger n, yields 

Cr-C = ^. (9) 

Therefore there is no multi affinity in this regime. The 
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FIG. 4: Ck{n,g) as a function of n for k = 1, 2 and 3. The 
three leftmost straight hnes have slopes according to Eq. 
while the bold middle line has a slope equal to 2/3 in ac- 
cordance with Eq. @. For large n the structure functions 
become flat indicating that one has reached the decorrelated 
regime. 

n for which there is the crossover between the Levy and 
the self-afRne regime will depend on k and is governed by 
prefactors that the scaling analysis presented here can- 
not access. For n beyond ^, the front decorrelates and the 
structure function becomes independent of n. We show 



in Fig. 01 the = 1, 2 and 3 structure functions. Their 
behavior is in accordance with our predictions. However, 
note that for k = 2, (2 = 7/12 = 0.58 which is close to 
C = 2/3. Furthermore, the self-affine regime is close to 
the decorrelated flat regime. Hence, it is hard to distin- 
guish between the Levy and the sclf-affinc regime for this 
value of k. As k increases, the Levy regime grows, as the 
overhangs are emphasized for larger k. 

To conclude, we have shown that the structure of the 
interface in a gradient percolation problem combines frac- 
tal and self-affine properties. The perimeter that includes 
numerous overhangs has the classical fractal structure 
[lo|. However, Solid-on-Solid fronts that arc extracted 
from the perimeter, have a clear sclf-affinc property up 
to a crossover length scale ^ even if local jumps, inherited 
from overhangs, are removed. On larger scales it shows 
an uncorrelated noise behavior. The structure function 
is, however, sensitive to the overhangs on smaller scales 
and this implies a multi-affine scaling behavior in this 
regime. Implications of our results for physical inter- 
pretations of analogical and numerical experiments are 
important. 
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